Finite doping signatures of the Mott transition in the two-dimensional Hubbard model 
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Experiments on layered materials call for a study of the influence of short-range spin correlations 
on the Mott transition. To this end, we solve the cluster dynamical mean-field equations for the 
Hubbard model on a plaquette with continuous-time quantum Monte Carlo. The normal state 
phase diagram as a function of temperature T, interaction strength U and filling n reveals that 
upon increasing n towards the insulator, there is a surface of first-order transition between two 
metals at non-zero doping. For T above the critical end line there is a maximum in scattering rate. 

PACS numbers: 71.27.+a,71.30.+h 



Simple models in physics can have surprisingly rich sets 
of solutions that explain complex behavior observed in 
nature with minimal assumptions. The Hubbard model 
is one such model. It contains a kinetic energy term that 
represents the band structure, and an on-site repulsion. 
That model has helped understand for example the Mott 
transition, that is the metal insulator transition induced 
by strong screened Coulomb interactions between elec- 
trons, a fundamental phenomenon in condensed matter 
physics [UH]. Several transition metal oxides, for exam- 
ple, are close to a Mott transition and their properties 
can be explained by the Hubbard model using single-site 
dynamical mean-field theory (DMFT) [3H5]. 

However, for layered doped Mott insulators, such as 
high temperature superconductors, or for layered organ- 
ics that show unusual Mott critical behavior [6], it is 
known that single-site DMFT fails because it does not 
take into account short-range spin correlations. Clus- 
ter generalizations of DMFT [U [7] are computationally 
expensive but they cure this shortcoming. Recent stud- 
ies have explored the interaction-driven Mott transition 
using cluster methods, indeed revealing sharp modifica- 
tions to the single-site DMFT picture [SJ [§]. However, 
an overall perspective on the transition induced by dop- 
ing is missing since previous work has focused on large 
interaction strength [TDlU5| . 

Here we report the whole normal-state phase diagram 
of the two-dimensional Hubbard model as a function 
of chemical potential /z, temperature T and interaction 
strength U . Our results are summarized in Fig(T] that 
we shall explain below. This phase diagram could also 
be explored with optical lattices of cold atoms where U, 
T and n can be controlled [THl HZ]. In short, for suf- 
ficiently large U, there is a first-order phase transition 
between two paramagnetic metallic phases, one of which 
evolves continuously from the Mott insulator. There is 
a large scattering rate associated with that transition. 
This maximum in scattering extends above the critical 
end line that terminates the first order surface at pro- 



gressively lower T and larger doping with increasing U . 
The broad maximum in the scattering rate is centered 
around what is optimal doping for the high-temperature 
superconductors. 

Model and method: The Hamiltonian for the two- 
dimensional Hubbard model reads 



H =~Y1 ^Aa^a + U X ( - 9 ) ( n if - 9 ) "A* X/' 



where Ci a and cf a operators annihilate and create elec- 
trons on site i with spin a, and n — cf a c ia . t is the 
hopping amplitude between nearest neighbors, U is the 
energy cost of double occupation at each site and fi is the 
chemical potential. 

We solve this model using cluster dynamical mean-field 
theory. The lattice problem ([I]) is mapped onto a 2 x 2 
plaquette immersed in a bath that is determined self- 
consistently in such a way that infinite lattice and pla- 
quette have the same self-energy. The partition function 
of the plaquette coupled to the bath is given by [4j 

Z = J V[l/>i,1p] e~ S "-fo dT/,f dr'E K ^(r)A K (r,r')V*(T') ) 

(2) 

where S c is the action of the cluster, K label the cluster 
momenta, and A is the bath hybridization matrix. The 
self-consistency condition reads 



A(iw n ) =iu n +fj,- ^ c (iuJ„) 



iu n + fi- t c (k) - S c (iw„) 



(3) 



Here S c is the cluster self-energy matrix, t c {k) is the 
matrix of hopping in the supercell notation and k runs 
over the reduced Brillouin zone of the supcrlatticc. 

To solve this problem, we use the numerically exact 
continuous time quantum Monte Carlo method [HI [19] , 
a finite temperature approach that relies on the Monte 
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Carlo summation of all diagrams obtained from the ex- 
pansion of the partition function ^ in powers of the 
hybridization A. This method does not have errors asso- 
ciated with time discretization or bath paramctrization. 
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FIG. 1: (a): Chemical potential-interaction-temperature 
phase diagram of the two-dimensional Hubbard model ob- 
tained by cluster DMFT. Cross-sections at constant U are 
shown. Dark-gray (blue) shaded regions represent the coex- 
istence of two phases. Light-gray (yellow) areas denote the 
onset of the Mott insulator state, characterized by a total oc- 
cupation n — 1. Projections on /i = 0, and T = planes are 
also shown (full lines and dashed lines respectively). Open 
dots mark the extrapolated T — values of U c i and U C 2- A 
critical end line (dotted line) begins at the full dot Umit- 
(b),(c),(d),(e): the cross-sections at constant U are shown in 
several T — \i plots. Three phases can be distinguished: metal 
(open circles), Mott insulator (crosses), momentum selective 
phase where some regions of the Brillouin zone are gapped 
and others gapless (filled squares). 

Phase diagram: Our results are summarized by the 
temperature T versus interaction strength U and chemi- 
cal potential fi phase diagram shown in Fig. [T^,. We per- 
formed scans at constant U, varying the chemical poten- 
tial for several temperatures. These scans are indicated 
as cross-sections in the phase diagram and as distinct 
T — [i planes in Figs. [TJs to[Tj3. Due to particle-hole sym- 
metry, the phase diagram is symmetric with respect to 



the [i = plane, where the model is half-filled, fn that 
plane, we find the results of Ref. [8] where the U driven 
first-order metal-insulator transition (MIT) is bounded 
by the spinodals U C \(T) and U C 2(T) (dashed lines) be- 
tween which two solutions of the cluster DMFT equa- 
tions, one metallic-like, and one insulating-like, coexist. 
The spinodals end at a second order point U mit ■ 

We found that the coexistence region extends to finite 
values of the chemical potential ji, and is delimited by 
the spinodal surfaces /j, c i(U, T) and /z c2 (£7, T) that end 
at a critical line (dotted line). The intersections be- 
tween these surfaces and the cross-sections at constant 
U are displayed in the phase diagram as dark grey (blue) 
shaded areas. For scans at U — 5.6t which lie in the 
region U cl (T = 0) < U < U c2 (T = 0), only the /Lt cl line 
exists. For scans at U — 5.8t corresponding the region 
U C 2{T = 0) < U < Umit, we find coexisting solutions at 
finite doping between fi c i(T) and fi C 2(T). As U increases, 
this region narrows and occurs at lower temperature. 

The Mott insulating phase, characterized by a plateau 
in the occupation at n — 1, is shown for various constant 
U planes as light grey (yellow) shaded regions. Even 
though in the region U cl {T = 0) < U < U c2 (T = 0) the 
Mott insulator can coexist with a metal, as indicated by 
a different shade of grey (yellow), this is generally not so. 
In all cases displayed, the /j, c2 (U,T) surface marks the 
vanishing of a metallic solution while the fj, c i(U,T) sur- 
face does not coincide with the vanishing of the Mott in- 
sulating state. Doping instead occurs gradually in certain 
cluster momenta, a phenomenon called 'sector-selective 
transition' |T3l [14]. This feature is the coarse grained 
analogue of a continuous appearance of the Fermi surface 
with doping, supporting the formation of arcs or pockets 
in the Fermi surface at small doping. 

More precisely, in the plaquette scheme there are four 
cluster momenta, (0,0), (tt, Q), (0,7r) and (7r,7r), which 
can be thought to represent a coarse grained approxi- 
mation of the Brillouin zone. In the region between the 
Mott insulator and \i c \ , the (tt, tt) sector remains gapped 
and its occupation does not change with /i, but the other 
sectors are gapless (cf filled squares in the T — fj, cuts of 
FigsflJ) to [TJj). For \i < ^ ci , all the sectors are gapless. 
The selective metallization of the Mott insulator is not 
restricted to moderate values of the coupling, but also 
takes place for large values of U. The transition between 
the sector selective metal and the metal can be first order 
or it can occur as a crossover depending on U, T and \i. 

First-order transition at finite doping: The plot of the 
occupation n as a function of /i, shown in Fig. [2] for dif- 
ferent values of U and a wide range of T is the best way 
to exhibit the first-order transition. A plateau in the 
curves at n(ji) — 1 appears above the critical coupling 
U c \ (T) and reveals the onset of the incompressible Mott 
state. The curves n(/z) display a strong temperature de- 
pendence, especially in the phase nearest to the Mott 
insulator, indicating that the strong correlation effects 
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FIG. 2: Occupation n versus [i for different values of in- 
teraction strength U. The data shown are for temperatures 
T/t = 1/10,1/25,1/50 (triangles, squares and circles respec- 
tively). When two solutions are found to coexist, the solutions 
obtained following the metallic and the insulating solution are 
indicated as full and open symbols respectively. The plateau 
in the occupation at n = 1 signals the onset of the incom- 
pressible Mott state. 



appear at low energies. 

For U/t = 5.6,5.8,6 and low temperatures, we find 
an hysteresis loop between two solutions, an unambigu- 
ous signature of the first order character of the transi- 
tion. The Clausius-Clapeyron equations reveal that the 
sector selective metal closest to the Mott insulator has 
smaller entropy and smaller double-occupation than the 
large doping phase. The hysteresis loop becomes larger 
with decreasing T. For U — 5.6t, located in the region 
U cl (T = 0) < U < U c2 (T = 0), the branch of the hys- 
teresis cycle that corresponds to increasing \x (i.e. the 
metallic branch), continues up to half-filling, at [i = 0. In 
contrast, for U larger than U C 2(T = 0), namely U = 5.8i 
and U = 6t, the metallic branch shows an upward jump 
at finite doping. The metallic branch endpoint defines 
fi C 2- On the other hand, the upper branch of the hys- 
teresis cycle, obtained by decreasing [i from the insula- 
tor, evolves continuously into the sector selective metal 
and then undergoes a downward jump at finite doping. 
This defines \i c \. 

The first order transition occurs between two different 
metals. Since these phases have the same symmetry, the 
first-order surface can end at a critical line at finite tem- 
perature. The n(n) curves for U — It do not show an 
hysteresis loop but the jump in compressibility (dn/dfijT 
at the lowest temperature, still visible, suggests that the 
hysteresis loop is very small and that we are near the crit- 
ical line. In addition, the crossing of the isotherms near 
n ~ 0.91 translates into a vanishing expansion coefficient 
(dn/dT)^ = 0. That feature, in this temperature range, 



continues to larger doping as U increases (not shown) as 
a left-over of the critical end line. 

Discussion: A natural way to understand the effects of 
the short-range spin correlations on the Mott transition, 
is to compare our cluster DMFT phase diagram, FigjTJi, 
that accounts for these magnetic correlations, with that 
of single-site DMFT, where those correlations are absent. 
In single-site DMFT, the half-filled Mott insulator is con- 
nected to a correlated metallic state through a first order 
transition [3J, which can take place as a function of either 
U, T or /i. The first order surface ends at finite temper- 
ature at a critical line . 

Short-range spin correlations modify the single-site 
picture. In the [/-driven MIT, scrutinized in Ref. [8], 
and visible in our phase diagram along the plane n = 0, 
the transition is still first order, but the critical value of 
the Mott endpoint is reduced by roughly a factor of two 
and the shape of the spinodal lines is modified because 
the T = insulating state has no entropy when magnetic 
correlations are taken into account. Our exhaustive cal- 
culations in the T — [i — U space reveal that the short 
range spin correlations have an additional surprising ef- 
fect. They shift the spinodal surface \i c \ to finite doping, 
i.e. the first order transition for U > U C 2(T) is between 
two metals instead of between an insulator and a metal. 

Insight on the role of spin correlations can also be ob- 
tained from the plaquette eigenstates [TT]. We find, at 
low temperature, that the states with largest probability 
are the singlet with 4 electrons in the cluster momentum 
K = (0,0), the half-filled triplet in K = (n, ir), and the 
doublet with 3 electrons and K = (tt, 0). In Figs. [3^ to 
[3ji we show their statistical weight P as a function of /x. 
At low doping, the half-filled singlet dominates, revealing 
that the spins are bound into singlets due to the superex- 
change mechanism. We also find that the hybridization 
function depends only weakly on the sector, consistent 
with spin correlations within the plaquette that domi- 
nate. At large doping (small /i), the charge excitation 
becomes more important, as expected, so the weight of 
the doublet at N — 3 has the highest probability. 

For large values of U (see U — 7t), at this T the weight 
of the doublet is larger than the triplet in the whole dop- 
ing range and no coexistence is found. On the contrary, 
for intermediate values of U (see U/t = 5.6,5.8,6), the 
doublet and the triplet states compete. The coexistence 
region is associated with this interplay: in the metallic 
branch of the hysteresis loop, the triplet prevails over the 
doublet; in contrast, in the insulating branch, it is the 
doublet that wins over the triplet (similarly to the large 
U case at low doping) . Since the doublet at N = 3 corre- 
sponds to charge excitations, and the triplet at N = 4 to 
spin excitations, the first order transition is associated to 
a competition between spin and charge. Once the spin- 
charge competition is removed, at sufficiently large U for 
our lowest T, we are above the critical end line and a 
smooth crossover takes place. 
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FIG. 3: (a) to (d): statistical weight P of the following 
cluster eigenstates as a function of n, for several values of 
U, and at the low temperature T/t = 1/50: the singlet 
\N = 4,5 = 0,K = (0,0)}, triplet \N = 4,5 = 1,K = {tt,tv)}, 
and doublet \N = 3, 5 = 1/2, K = (7r,0)} (squares, triangles 
and circles respectively) where N, S and K are the num- 
ber of electrons, the total spin and the cluster momentum 
of the cluster eigenstate. (e) to (h): scattering rate r=- 
ImE(^ ! o)( tJ — ► 0) in the sector (n,0) versus /x, at T/t = 1/50. 



U = I2t in Ref. [TTJ The crossover, as measured by 
the maximum scattering rate or by (dn/dT)^ = 0, is 
caused by the influence of the critical line that moves 
to progressively lower temperature and higher doping as 
U increases. Accounting for differences in methods, this 
crossover caused by critical behavior is also the one de- 
tected for U = St in Refs. UHH We cannot strictly 
exclude that the critical line that ends our first-order 
surface becomes a quantum critical line or point at some 
critical U but this possibility is unlikely. 

We conjecture that the quantum critical behavior sur- 
mised in high-temperature superconductors [231 124] is a 
constant U cut of the low temperature behavior identified 
in our work as originating from the finite-doping critical 
end line of the first-order surface between two types of 
metals. That transition occurs only when U is larger 
than the critical value U c i for the Mott transition. The 
influence of Mott physics extends way beyond half-filling. 

We are indebted to S. Allen for technical help and to D. 
Senechal for useful discussions. This work was partially 
supported by FQRNT and by the Tier I Canada Research 
Chair Program (A.-M.S.T.). Computational resources 
were provided by CFI, MELS, the RQCHP and Compute 
Canada. 



Competing charge and spins fluctuations near the Mott 
transition are a source of scattering, as can be seen 
from Figs. [3^ to [3ja. There we plot the scattering rate 
T{n,T/t = 1/50), estimated from the extrapolation of 
the imaginary part of the cluster self-energy in the (tt, 0) 
sector. In a Fermi liquid, Y as well as the scattering rate 
at all other cluster momenta, should vanish. In contrast, 
leaving aside the insulator, we find a maximum in T ei- 
ther near the hysteresis loop of the first order transition, 
or near the critical end line at large U (see for example 
U = It). We have verified that as temperature increases, 
the value of T at its maximum increases as docs its width 
in doping, reminiscent of the behavior observed around 
optimal doping in the cuprates. 

Several groups, using different cluster methods, al- 
ready reported an anomalous metallic state at finite dop- 
ing, linking this to the pseudogap phenomenon [TU1 HH 
ITS] , to the selective transition [T31 E] , or to a competi- 
tion between Kondo and superexchange scale [TT1 ITU I2T] . 
However, those investigations mostly focused on a region 
where the interaction strength is large. Our contribution 
is to track the origin of the large scattering rate to the 
influence of the Mott phenomenon at finite doping, where 
spin and charge fluctuations compete. 

Although the first order transition and its associated 
maximum in Y occurs at small doping in the T > O.OTi 
range that we can access, its signature in the form of 
a crossover occurs at progressively larger doping as U 
increases, reaching the optimal doping point found for 
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